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The average number of solutions 
of the Diophantine equation = W'^ 

and related arithmetic functions 

M. Kiihleitner and W.G. Nowak 



Abstract. For the number of integer solutions of the title equation, with W < x {x a. large parameter), 
an asymptotics of the form Ax\ogx + Bx + O ^a;^/^(log a;)^(log log a;)^) is established. This is achieved 
^ I ' . in a general setting which furnishes applications to some other natural arithmetic functions. 

^ ! AMS-Classification. 11N37, 11M06, 11R42, 11D25 

■3 

1. Introduction. During tlie problem session of the 1991 Czechoslovak Number 
Theory Conference, A. Schinzel proposed the following question: Let r(n) denote the 
number of ways to write the positive integer n as a sum of two squares, and consider 
the asymptotic formula 



^{r{n)f = 4:xlogx + Cx + o(x^/^u{x)^ , (1.1) 



as X oo. How small can the factor uj{x) be made, if one uses the sharpest tools of 
^ ■ contemporary analytic number theory? 

j> I In fact, asymptotics of the shape (1-1) had been established by B.M. Wilson [19], 
with uj{x) = x^ , and W. Recknagel [13], with uj{x) = (logx)^. 



\-i \ Subsequently, as a reaction to Schinzel's proposal, the first named author [6] sharpened 
' the estimate to uj{x) = (logx)^^/'^(loglogx)^/'^. 

The problem addressed by Schinzel is closely related to another arithmetic question to 
which our present title refers: Given a large parameter x, how many integer triples 
(w, f , w) exist with u'^ +v'^ = , w < x 7 This matter was dealt with by K.H. Fischer 
[2] and also by Recknagel [13]. 

Since ^r(-) is multiplicative, it is easy to write up the corresponding generating Dirichlet 
series for (r(?T,))^ and r{n^) . (See section 4.2 for details.) Thus, both problems are 
subsumed in a natural way by the following more general result which is to be the 
objective of the present note. 



^^"^ To complete the history, one should mention the older, somewhat coarser results of W. Sierpihski 
[16] and S. Ramanujan [12], with error terms O [x^/'^logx) , resp., Furthermore, Schinzel 

[15] himself had bounded the remainder from below, showing that it is an n{x^/^) . 



2 



Theorem. Let a{n) be an arithmetic function satisfying a{n) <^ for every s > , 
with a Dirichlet series 

where t/ie Dedekind zeta-function of some quadratic number field JK , G{s) is 

holomorphic and bounded in some half-plane 3?(s) > 9 , 9 < ^ , and mi , 777-2 ore 
nonnegative integers. Then, for x large, 



n<x 



J2 <ri) = Res (^F{s) y j + O (x^/^{\ogxf{\og\ogx)"''+"'^^ = 
= Axlogx + Bx + O (x^/^{\ogxf{\og\ogx)'^'+"'^^ . 



Remark. It is natural to ask why we have chosen just this very degree of generality, 
resp., speciahzation, in our suppositions. Note that, for IK a quadratic field, Ck(-s) = 
C{s)L{s,x) 1 where x is a certain real Dirichlet character. Thus the above function F 
is a special case of 

L(.,xW)...L(s,x(^^) . . 
Ck,(2.)...Ck^(2.) 

where the L{s,x^''^) ^re ordinary Dirichlet L -functions, IKi,...,IKm are arbitrary 
algebraic number fields, and G{s) is as before. As long as J = 4, our argument 
would go through without major alterations and lead to the same result, apart from the 
exponents of the log- and loglog-factors in the error term. 

However, for J 7^ 4 , the situation changes drastically: If J = 2 or 3 , an important role 
is played by the zero-free region of the denominator. Using the sharpest information 
available of this kind, one obtains a bound 0(a;^/^ exp(—c(loga;)"^/^(log log for 
the remainder, where the exp-factor is familiar from the prime number theorem. (See, 
for instance, formulas (14.28) - (14.30) and (1.105) - (1.107) in Ivic [5].) 
Opposed to this, for J > 4 the size of the nominator in vertical strips becomes 
important. This ultimately leads to an estimate 0{x'^-^~^'^) , where aj is the best known 
error-exponent in the J-dimensional Piltz divisor problem which, for J > 5, is still a 
matter of small improvements from time to time (cf. Ivic [5], p. 355, or Titchmarsh [18], 
ch. 12). Thus the case J ~ A remains as the most delicate one, at least as far as smaller 
factors, apart from x^/^ , in the error term are concerned. 

In addition, our analysis has been inspired by recent work^^^ of K. Ramachandra and 
A. Sankaranarayanan [11] which deals with 

^^•^ Therefore, in these cases, the assumption of the truth of the Riemann Hypothesis (RH) leads to 
a better error term 0{x^) , with some 6 < ^ . In contrast, our Theorem is as sharp as it would be if 

RH could be proven. 

The authors are indebted to Professor Sankaranarayanan for sending them a preprint of this paper. 
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n=l 



See also Sankaranarayanan's Oberwolfach lecture [14]. However, our argument is tech- 
nically a bit simpler. 



2. Some auxiliary results. 

Lemma. Let IK be an arbitrary number field of degree [JK : Q] > 1 , its Dedekind 

zeta- function, and e > fixed. Then, for each sufficiently large T there exists a 
measurable set At — At,k C [T, 2T] of Lebesgue measure 

X{At)<T', (2.1) 

with the property that, for some Cs^k > depending only on e and IK , 

sup |Ck(1 + ^^)^' < Cs,K loglogT. (2.2) 

te[T,2T]\AT 



Proof. For the case of the Riemann zeta-function, this is contained in Ramachandra 
[10], Theorem 1. Our argument follows the lines of Lemma 3.2 in Ramachandra and 
Sankaranarayanan [11]. 

According to Heath-Brown [3], there exists some 5 = 5{e) > such that, for T large, 
the cardinality of the set 

7^,,T = {peC: CM = 0, 5ft(p) > 1 - 35, T < Q{p) < 2T } (2.3) 

satisfies 

#7^,,T < T'/^ . (2.4) 

For any fixed c > 0, we define 

^(c,T) := [T,T + cT^/4]U[2T-cT^/^2T]U |J [^(p) - cT^/^ $>(p) + cT^/^] , 

and 

M{c,T) := {s e C : U{s) > 1-35, ^{s) e [T, 2T] \A{c, T) }U{s e C : ^{s) > 1 } . (2.5) 

On every M.{c,T), by construction Ck(s) ^ 0, thus logC]K;(s) can be defined properly 
by analytic continuation . For any r] = 2 + it with t e [T, 2T] \A{1,T), we consider 
the circular discs 

Ci := {s e C : \s-r]\<l + 5}, C2 := {s E C : \s - r]\ < 1 + 25 } . 

^^■^ To be explicite, on A4{c,T), we can define log(;'][^(s) := log(^]g(2) + _/^(Ck(-^)/Ck(^)) ^-^ where C 
consists of the two straight Hne segments from 2 to 2 + G(s)i and from 2 + 'ts{s)i to s . 
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Evidently, Ci C C2 C M{^,T). Thus the Borel-Caratheodory inequahty (cf. Titch- 
marsh [17], ch. 5.5) yields 

max|logCK(s)| < T maxlog|CK(s)| + — ^ — |log Ck(2 + . 

seCi s€C2 

C 1 

Since, for a certain C > 0, |Ck(o" + ^^)| < |^| , provided that a > ^, \t\ sufficiently 
large, it follows that log|CK(-s)| < C'log(2T) on C2. Hence, in particular, 

\\ogCK{(T + it)\<^logT for (7 > 1-5, 

the constant involved not depending on t. Therefore, 

logCK(s)<logT (2.6) 

uniformly for all s with ^{s) >1-S, '^{s) G [T, 2T] \ ^(1, T) . 
We now use the series representation (valid for 3?(s) > 1 ) 

CO 

iogCK(«) = E E = T.^i'pr' + H{s) , (2.7) 

1 ft 

V m=l V 

where H[s) is regular and bounded in any half-plane ^{s) >e>\. {V are the prime 
ideals in the ring Ok of algebraic integers in the field IK and A/'(-) denotes the norm 
of ideals in O-^.) 

We put X = (logT)V'5. Then, for any t e R, 



(5+icxD 

1 



27ri 

5 — ioo 



log Ck(1 + it + w)T{w)X'^ dw 



00 I '^+*°° \ 

E E ^ / (^"'^(^)'")""r(«;)d«; 



CO ^ 

= ^ 5^ _AA(p)-^(i+^*) exp(-A/'(7^)"^/X) = 5^Ar(7')-^-'*exp(-A/'(P)/X)+0(l) . 

•p m=l V 

(2.8) 

To estimate the remaining sum, we observe that 

^(p)-i-*exp(-A^(7^)/X)« M{^)-^^M-m)/X) = 

N{V)>X Ar{J)>X 



_^ \A/'(X)<u 



(5) 



Note that, by (2.7), |log Ck(2 + < logCK(2)- 



using only that ^ 1 <^ u. (I denotes arbitrary integral ideals in Ck-) Further, 
AA(X)<u 



J2 M{V)-^-''eM-^f{V)/X) 



< ^(^)~' = loglogX + 0(l), 

U{V)<X 



where the last conclusion is immediate from the prime ideal theorem in the form 
(cf. Narkiewicz [8], pp. 369-372) 



E 1 = 



u 



+ 



log u \ (log u) 



u 



N{V)<u 

In view of (2.8) and the choice X = (logT)-*^/^, this implies that 



1 



(5+icxD 



— / logCK{l + it + w)T{w)X'^dw 



5—ioo 



< log log log r + 0(1). (2.9) 



Our next step is to put W = (log log T)^ and to recall Stirling's formula in the weak 
form (valid uniformly in any strip cri<cr<a"2, |t|>l) 

T{a + it) < |t|''"^/^exp(-f |t|). (2.10) 
From this it readily follows that 

(5±ioo oo 

J logCKil + it + w)T{w)X'^dw<^X^ J exp(-fu)dw<X'^exp(-fW^) < 1, 

5±iW 

hence 



w 



5+iW 



/ logCKil + it + w)T{w)X'"dw 
2-111 J 



S-iW 



<logloglogr + 0(l). (2.11) 



From now on we impose the condition that t G [T, 2T]\A{2, T) . We evaluate the integral 
in (2.11) by the residue theorem, applied to the rectangle K with vertices ±6 ± iW . 
(Note that for t G [T, 2T] \ A{2, T) and w G 3?, necessarily I + it + w e M{1, T) , and 
for log(CK(l + it + w)) the bound (2.6) applies.) We obtain 

-5+iW W 

j \ogCK{l + it + w)T{w)X'" dw ^X-^ XogT J |r(-5 + m)| dw < 1 , 
-S-iW -w 

and similarly, for the horizontal segments, 

5±iW 

J logCK{l + it + w)T{w)X'^dw<^X^ logTexp(-f^y) < 1. 

-S±iW 
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Since the only pole inside the rectangle, at w = 0, gives a residue of logCK(l + it) , we 
altogether derive from (2.11) that 

|log Ck(1 + it) I < log log log T + Ci , 

for all t e [T, 2T] \A{2,T) (Ci > an appropriate constant depending on K and s). 
Hence, for the same t, 

\Ck{1 + it)\~^ = exp(-log|CK(l < exp(|logC]K(l + < loglogT. 

Taking At,k = A{2,T), this just proves clause (2.2) of the Lemma, while (2.1) is 
immediate by (2.4) and (2.5). [| 

We conclude this section by stating some more bounds for the zeta-functions involved, 
which will be needed in the proof of the Theorem. First of all, again for every number 
field IK with [K : Q] > 1, 

|CK(a + zt)r'«(log(2 + |t|))^ (2.12) 

uniformly (at least) in a > 1 , the constant C > possibly depending on IK . (This is 
most conveniently deduced after the classic example of Apostol [1], Th. 13.7, using the 
necessary facts about the Dedekind zeta-function from Narkiewicz [8], ch. 7.2.) 

Further, we recall that 

C(| + it) « (1 + |t|)V6+- , L{-e' + it, x) «x (1 + 1*1)'/'+^ 

for any Dirichlet L-series. (For the first estimate, see Titchmarsh [18], Th. 5.12. The 
second one follows readily from the functional equation in Apostol [1], Th. 12.11, along 
with Stirling's formula.) Applying the Phragmen-Lindelof theorem in the form given 
by Titchmarsh [17], ch. 5.65, we infer that 

C(a + it)« 1 + 1^1^'-"^/'+^ , L{a + it,x)<^x l + |t|^'-"^/'+' , 

uniformly (at least) in cr > i, |t| > 1. Now, if IK is a quadratic field of discriminant 
D, its zeta-function can be factorized as Ck{s) = C{^)L{s,xd) ■ (Cf. Zagier [21], p. 100, 
in particular for the definition of the real character xd involved.) Therefore, in this 
case, 

CK(a + it)«l + |t|^('-")+' , (2.13) 

uniformly in o" > | , |t| > 1. Last but not least, it has been proved by W. Miiller [7] 
(even in the stronger form of an asymptotics) that 

T 

J |CK(i+it)|'dt«T(logT)2, (2.14) 



In fact, the estimates (2.12) and (2.13) are far away from being the sharpest ones of their kind, 
but they suffice for our purpose and are conveniently available in textbooks. 
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again for every quadratic field K . 



3. Proof of the Theorem. W.l.o.g., let x be half an odd integer. Then, by a simple 
version of Perron's truncated formula (see, e.g., Prachar [9], p. 376), 

for any fixed e > sufficiently small. We shift the path of integration to the vertical 
line 3?(s) = | . The horizontal segments contribute 

l+e±ix3/5 1+e o/ro2 



By (2.12) and (2.13), this is 

max (x^x'/'f^'-^^/')<^x^/'+' 

|<CT<H-e V / 

Hence, by the residue theorem. 



F(s)x'' — -C X 
s 



^ a(n) = Res (f{s) y ) + 3^ / ^(s) x'^ + O (x^/^^ . (3.2) 

To estimate the remaining integral, we use the parametrization s = ^{1 + it) , \t\ < 
2x^^^ , along with a dyadic decomposition, to obtain 

ds ^ 1 / |Ck(Ki + ^-^))| 

_4^3/. IC(l + ^^)r ICK(l + ^^)^ 

(3.3) 

where the sum in fact contains only O(logx) terms. With ^k,t as in the Lemma, we 
put i3K,T = [T,2T] \^K,T. By (2.1), (2.12), and (2.13), 

I KK(^(i + ^t))|^ dt-r^/^+^- 

Moreover, by (2.2) and (2.14), 

|Cic(i(l + ^t))|' 



/ IC(i 



+ ^^)r lCK(l + ^^)^^ 



dt<T(logT)2(loglogT) 



mi+m2 
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Inserting the last two bounds into (3.3) and summing over T, we obtain 

\mi+m2 



/ds 
F(s)a;* — -C x^/^(logx)^(loglogx)' 



Together with (3.2), this completes the proof of our Theorem. 

4. Applications. It remains to verify that the result established covers the two special 
problems addressed in the title and in the introduction, even in a more general context. 

4.1. The second moment of quadratic Dedekind-zeta coefficients. For an 

arbitrary quadratic number field IK with discriminant D, let Or denote the ring of 
algebraic integers in IK , and rK(n) the number of integral ideals X in Ok of Norm 
N'i'I) = n . Then, as we shall show below, for 3?(s) > 1 , 

f^i^j^Yli^^p-r- (4.1) 

n=l ^ p\D 

Therefore, applying the Theorem with (mi, 777,2) = (1, 0) , we obtain what follows . 
Corollary 1. For any quadratic field IK of discriminant D, and x large, 
^{f'^{n)f = ^ia;loga; + 5ia; + (^x^/^ (log a;)^ log log a;) , 



n<x 

with 



p\D 



p\D \ p\D ^ 

where 7 is the Euler-Mascheroni constant. 

Since r{n) = 4rQ(i)(?7) , this implies in particular that 

^(r(n))^ = 4a; log a; + IGBix + O ^a;^/^ (log a;)^ log log a; j (4.2) 

n<x 

with 

4 32 48 

I6S1 = -4 + 87 + - log 2 + —L'{1) - -^C'{2) ^ 8.0665 , 

3 TT TT 

where L(-) is the L -function corresponding to the non-principal Dirichlet character 
mod 4. 



^^•^ For the evaluation of the residue at s = 1 and the subsequent numerical computations, we have 
employed Mathematica [20]. 
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It remains to verify (4.1). To this end, we recall the decomposition laws in (cf. Zagier 
[21]). On the set IP of all rational primes p , there exists a partition IP = IP0UP1UIP2 
such that 

(i) peWo ^ p\D ^ {p)^V\ Af{V)^p, 

(ii) peiPi ^ (p) = riV2, Vi^V2, ^f{Vl)=^^{V2) = p, 

{in) p e P 2 ^^=^ {p) prime in O-^ , N'{p) — p^ . 

For 3?(s) > 1 , this implies the Euler product representation 

pep \ j=i ) V (4.3) 

= n (i-p-^)"' n (1-^^-^)"' n a-^^-'r^ 

peiPo peiPi peiP2 



By comparison, we immediately see that r^{jP) — 1 throughout for p G Pq, ?']k(p') 

1 if 2 is even, 
if j is odd. 



1 if 7 IS GVGn 

J ■ + 1 for p e Pi , and r^{p^) = \ ■ , , ' for p e P2 



Hence, 



E^^= 11(1 -f-r' n (i + Eo+DV^'l Ild-f-")-'- 

n=i peiPo peiPi \ j=i / p€]P2 

Using that ^^(j + \)'^z^ = (1 + -2)(1 — (for < 1), we see that this equals 
i>o 

n (i-p-^)-' n (i-p-^)-'(i+p-^) n (1-^-'^)-'= 

peiPo peiPi peiP2 

: n(i -^'"'^) ( n (1 -v-)-^ n ^^-v-^'^ n ) n {^^v-r^ 

per \peiPo peJPi peiP2 / pePo 



C(2s) 



which proves (4.1). 



Note that, on the basis of this partition, the Dirichlet character xd can be defined as follows: 
X(p) = if p e IPo , X(P) = 1 if P e IPi , and x(p) = -1 if p G IP2 • 



10 



4.2. Diophantine equations as addressed in the title. Again for K a quadratic 
field, we shall show that 



n=l 



where G{s) is holomorphic and bounded in every half-plane 3?(s) > ^ > | . Taking 
this for granted, we apply our Theorem with (mi, 777,2) = (1, 1) and derive the following 
consequence. 

Corollary 2. For any quadratic field IK of discriminant D , and x large, 



r^{n^) — A2x\ogx + B2X + O ^x^^^ (log a;) ^ (log log 

n<x 

with 

, 36L(1,xd)'G(1) 



X? 



36L(1,xd)^G(1) ^ 2L-(1,xd) .^o^^'^^) 2L'(2,xd) 24C'(2) 
-^2 = ttt:; ^ — —rr, ^ — 1 + 27 + 



n^L{2,XD) \L{1,xd) ' G{1) L{2,xd) tt^ 

where G{s) is given in form of an Euler product in (4.6) below . 

To apply this result to Diophantine equations, let Q = Q{u, v) = au^ + huv + cf ^ be an 
integral, primitive, positive definite binary quadratic form of class number l.*^^^) Then 
it is well-known (cf. Zagier [21], §§ 8 and 11) that 

rQ{n) -.^ #{{u,v) e : Q{u,v) = n } ^ u;d rK{n) , (4.5) 
re for D = -3, 

where IK = Q(v^), and^^^) a;^ = < 4 for i:> = -4, 

I 2 else. 

We thus can infer the following conclusion. 

Corollary 3. For every integral, primitive, positive definite binary quadratic form Q 
of class number 1, with discriminant D , and for large real x , 

^{(•u, V, w) El? : Q{u, v) = to^, w < X } = 
= 1+ E rQ{n^) = l+u;D Yl = 

l<n<.x l<n<x 

= uj d{A2X log X + B2X) +0 ^x^/^ (log x)^ (log log x)^^ , 

where A2 , B2 are as in Corollary 2 . 

^'^^ Of course, G'(l)/G(l) is most easily evaluated as -g^logG(s)| _^ , which transforms the products 
into sums of derivatives of logarithms. 

^^^^ Equivalently, the discriminant D = — 4ac is one of the ^' numeri idonei" {—3, —4, —7, —8, —11, 

-12, -16, -19, -27, -28, -43, -67, -163}. See, e.g., Hlawka/Schoiengeier [4], p. 92. 
^^^^ u!D is the number of automorphisms of the form Q or, in other terms, the number of units in 
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In particular, for the Diophantine equation of the title, we obtain 

i^{{u,v,w)eZ^:u'^+v^=w^,w<x}^l+ ^ r(n^) = 

l<n<x 

= 4:A2xlogx + AB2X + O ^a;^/^(loga;)^(loglogx)^j , 

with 

4^2 = 5^ ^ 0.9091 , 

= £gai f «^'('. - 1 + 2, 1 - ^^''^;.^-^) _ g^m) . 2.1715 , 

TT^ C \ TT G(l) C TT^ / 

where C = L(2,x-4) = ^(~l)''(2j + is Catalan's constant, and, by (4.6) below, 

i>o 

n n , 

p=lmod4 '^^ ' p=3mod4^ 

^ = --log2-6 y , ,-4 y -0.35876. 



It remains to verify (4.4). By what we noted earlier (right after (4.3)), for 3fJ(s) > 1, 



-js I _ 



Using that ^^(3^ + l)z^ = (1 + 22;)(1 — z)~'^ (for \z\ < 1), we may write this as 



F(s)= n (i-f-)-' n n (^-p-^-r 

pePo pePi ^ ^ ^ peiP2 



Recalling the Euler product representation of Ck{s) and Ck(2s) , as given in (4.3), along 
with that of C(2s) , we obtain after a brief calculation 

n (i^^ n MP-) na-^'-'T', (4.6) 

^^jp^ p ) ^^^^ ^^^^ 

where h{z) := (1 + 2z){l - zY{l - z'^)'^ = 1 + 0{\z\^) around z = 0. Thus G(s) 
possesses a Dirichlet series absolutely convergent for 3?(s) > |. This proves (4.4). 
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